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Shallow-fluid  equation  models  have  been  applied  to  several  Army  problem  areas, 
especially  those  associated  with  atmospheric  transport  and  diffusion.  The 
objective  of  this  report  is  to  present  the  development  of  a system  of  equations 
which  govern  the  behavior  of  material  atmospheric  layers  and  which  retain  the 
advantages  of  the  shallow-fluid  equations  but  are  more  generally  suitable  for 
atmospheric  analysis.  The  report  consists  of  three  major  sections.  The  first 
section  presents  a formal  derivation  of  what  therein  is  defined  as  the  "Materia! 
Layer  Equations."  The  second  section  addresses  some  of  the  practical  \ 
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0.  ABSTRACT  (cont) 

considerations  involved  in  applications  of  the  material  layer  equations. 
Although  the  presentation  of  numerical  algorithms  for  specific  applications  of 
the  material  layer  equations  is  beyond  the  scope  of  this  report,  the  third 
section  does  present  elementary  examples  of  how  the  material  layer  equations 
may  be  specialized  to  address  specific  atmospheric  analysis  problems.  It  is 
concluded  that  the  material  layer  approach  is  most  useful  when  the  matter  of 
concern  is  the  transport  and  diffusion  of  material  such  as  from  NBC  weapons  or 
smoke. 
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INTRODUCTION 

The  Army's  requirements  for  meteorological  information  are  diverse  [1]. 
Over  the  last  several  decades,  the  foremost  Army  meteorological  problem 
areas  have  been  tube  artillery  ballistics  and  chemical /biological 
warfare/defense.  The  latter  is  especially  dependent  on  neteorological 
information  since  atmospheric  transport/dif fusion  is  an  integral  part 
of  weapon  deployment  systems  (similar  remarks  pertain  to  the  deployment 
of  smoke  as  a battlefield  obscurant). 

The  requirements  for  meteorological  information  cannot  always  be 
satisfied  by  one  or  even  several  atmospheric  observations.  The  atmo- 
sphere may  respond  to  complex  terrain  in  unexpected  ways.  This  problem 
was  manifested  by  the  alleged  "sheep"  incident  at  Dugway  Proving  Ground 
a decade  ago.  This  situation  motivated  investigations  regarding  the 
numerical  modeling  of  atmospheric  flow  in  association  with  complex 
terrain.  The  primary  approach  that  has  been  used  is  based  upon  the 
shallow-fluid  equations  (see  Meyers  and  Dettling  [2],  Tingle  and 
Bjorklund  [3,4],  and  Dumbauld  and  Bjorklund  C 5 ] ) . The  shallow-fluid 
equation  models  have  reached  the  stage  of  development  where  they  have 
had  operational  use  in  Army  test  and  evaluation  activities  associated 
with  transport  and  diffusion  problems.  On  yet  another  front,  a shallow- 
fluid  equation  model  has  been  developed  [6]  with  the  object  of  short- 
term mesoscale  predictions  for  artillery  ballistics. 

As  useful  as  the  shallow-fluid  equations  might  seem  for  some  Army  re- 
quirements, they  may  be  better  suited  to  the  Navy  since  the  equations 
are  far  more  valid  for  the  hydrosphere  than  the  atmosphere.  Acknowledg- 
ing the  utility  of  the  shallow-fluid  equations  but  recognizing  their 
limitations,  this  report  develops  the  theoretical  basis  for  "atmospheric 
material  layer  equations"  and  discusses  the  more  general  utility  of  this 
approach  to  atmospheric  problems. 


The  report  consists  of  three  major  sections.  The  first  section  presents 
a formal  derivation  of  what  herein  is  defined  as  the  "Material  Layer 
Equations."  The  second  section  addresses  some  of  the  practical  consid- 
erations involved  in  applications  of  the  material  layer  equations. 

Although  the  presentation  of  numerical  algorithms  for  specific  applica- 
tions of  the  material  layer  equations  is  beyond  the  scope  of  this 
report,  the  third  section  does  present  elementary  examples  of  how  the 
material  layer  equations  may  be  specialized  to  address  specific  atmo- 
spheric analysis  problems. 

MATERIAL  LAYERS 

The  atmosphere  is  a relatively  thin  layer  of  gas  clinging  to  the  earth 
as  a consequence  of  the  forces  of  gravity  and  friction  at  the  earth's 
surface.  The  atmospheric  layer  is  itself  internally  laminate,  con-  , 

slsting  of  overlapping  layers  with  differing  air  mass  properties. 

Observers  can  perceive  material  layers  in  the  atmosphere  when  ascending 
or  descending  in  an  aircraft  through  haze  layers.  These  layers  are 


especially  manifest  in  industrialized  regions  where  atmospheric  pollu- 
on  8 trapped  within  a material  layer  which  is  separated  from  a 
->uper  or  unpolluted  air  mass  by  a stable  subsidence  inversion. 

loimally,  a material  layer  is  one  which  is  mass  conservative  as  a con- 
sequence of  being  bounded  by  material  surfaces  through  which  there  is 
no  mass  transport.  A material  (or  Lagrangian)  surface  (interface)  is  one 
wherein  air  parcels  are  confined  to  two-dimensional  motion.  The  atmosphere 
does  not  contain  material  interfaces  such  as  those  which  occur  between 
immiscible  tluid  layers;  however,  there  are  hydrostatically  stable  layers 
in  the  atmosphere  which  commonly  behave  as  the  equivalent  of  immiscible 
material  layers. 


The  behavior  of  the  atmosphere  can  be  fully  understood  only  if  the 

rl!n,rt^nCe  ?f  maC®rlal  layers  is  adequately  recognized.  The  kinematic, 
1nd  waather  structure  of  the  atmosphere  is  intimately 
reiated  to  the  dynamics  of  material  layers.  One  of  the  greatest  break- 
throughs in  meteorology  has  been  the  Introduction  of  frontal  and  air 
mass  anaiysis  by  J.  Bjerknes  and  collaborators  [7],  leather  satellite 
pictures  continue  to  confirm  the  great  foresight  of  these  Norwegian 
Pi  neers  in  describing  the  properties  of  air  masses  (material  layers) 
and  associated  fronts  (material  interfaces) , and  their  role  in  extra- 
tropical  cyclone  development.  Frontal  inversions  are  only  one  example 
material  interfaces  which  affect  atmospheric  behavior.  The  tropo- 
pause,  subsidence  inversions,  and  the  upper  (nonturbulent)  portion  of 

ZlltCe  irT  arC  at,dltlonal  significant  atmospheric  material  sur- 
foremo8timaterlal  surface,  however,  is  the  air/earth  inter- 
face which  causes  major  terrain  and  diabatic  effects  on  atmospheric 
structure  and  behavior. 


8 *he  last  t”°  deca<*es,  a creeping”  breakthrough  has  occurred  in 
eteorology  associated  with  the  development  and  improvement  of  high- 
speed digital  computers  and  the  concomitant  development  of  numerical 

Propreoa°y  df?fn' ^stic/prognostic  analysis  of  atmospheric  structure. 
Progress  in  this  area  has  been  rapid  and  the  operational  primitive- 
equation  model  described  by  Shuman  and  Hovermale  I 8 1 is  surely  a great 
achievement.  Further  progress  is  assured;  the  only  practical  limit  is 
the  capacity  and  capabilities  of  computers. 

Unwittingly , the  numerical  model  fad  has,  especially  in  the  early 

.t0  rh  interest  in  material  layers  into  the  background. 
The  bold  atmospheric  structure  associated  with  material  layers  was 
incongruous  with  demands  for  smoothness  by  the  elementary  numerical 
models.  Fortunately,  there  is  increasing  relaxation  of  this  stricture- 
for  example,  the  primitive-equation  model  of  the  National  Weather 
pamtioningreCORnlZeS  malntainS  the  troposphere/stratosphere 


In  another  area,  in  the  last  several  decades,  interest  has  developed 
concerning  the  existence  and  nature  of  hydraulic  jumps  in  the 
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atmosphere  associated  with  the  flow  of  stratified  (material)  layers. 

On  the  one  hand.  Freeman  (9|  and  Tepper  ( 10 | have  suggested  that 
hydraulic  Jumps  associated  with  traveling  gravity  waves  on  subsidence 
inversions  are  the  long  sought  explanation  for  squall  lines.  On  the 
other  hand,  Houghton  and  Kasahara  111  I and  Long  I 12 1 have  mathemat icallv 
investigated  the  properties  of  hydraulic  jumps  associated  with  strati- 
fied flow  over  barriers. 

The  shallow-fluid  equations,  prominent  in  hydrodynamics  and  oceanography 
(see  Stoker  l 13  1 ) , are  the  approach  frequently  used  in  the  investigation 
of  the  flow  of  stratified  fluids.  Stoker  (14  | , kasahara,  Isaacson,  and 
Stoker  | 15  1,  and  Turkel  I 16  I have  used  various  forms  of  the  shallow- 
fluid  equations  to  investigate  frontal  motion  in  the  atmosphere. 

Freeman  [6  1,  Lavoie  { 1 7 1 , and  Tingle  and  Bjorklund  I 3,4  I have  reported 
on  numerical  shallow-fluid  models  which  are  useful  for  investigating 
terrain  effects  on  atmospheric  flow. 

Hydraulic  jumps,  if  only  transitory,  are  a characteristic  feature  of  the 
shallow-fluid  equations.  The  use  of  the  shallow-fluid  equations  in  the 
development  of  numerical  models  has  required  intensive  investigations 
of  numerical  algorithms  which  are  computationally  stable  and  yet 
accommodate  the  jump  discontinuities.  The  works  of  Lax  and  Wendorf  l 1 8 1 
and  Sundstrom  ( 19  I are  significant  examples  of  this  effort.  Although 
the  subject  of  hydraulic  jumps  and  their  solution  has  developed  into  a 
minor  fad,  this  does  not  detract  from  the  fact  that  the  introduction  of 
the  shallow-fluid  equations  into  atmospheric  dynamics  has  been  a major 
step  towards  returning  material  layers  to  the  prominence  that  they 
deserve  in  atmospheric  diagnostic/prognostic  analysis. 

The  classical  shallow-fluid  equations  pertain  to  the  flow  of  immiscible 
layers  of  incompressible  fluid  that  are  so  shallow  that  the  vertical 
shear  structure  is  essentially  irrelevant  and  the  mean  velocity  in  the 
layer  is  sufficient  to  describe  the  flow.  The  great  power  of  the 
shallow-fluid  equations  is  manifest  in  their  strong  conservation  of 
mass  and  energy;  however,  the  classical  shallow-fluid  equations  are 
only  crude  approximations  when  applied  to  the  atmosphere. 

MATERIAL  LAYER  DYNAMICS 

The  development  of  the  material  layer  equations  begins  with  a standard 
set  of  atmospheric  variables.  Initially,  the  coordinate  system  is  a 
Cartesian  conformal  projection  of  a relatively  small  region  of  a 
rotating  earth,  but  a vertical  coordinate  transformation  conforming  to 
the  material  surfaces  alters  the  basic  dynamic  equations  so  as  to  be 
consistent  with  the  hypothesis  of  no  mass  transport  through  the 
bounding  material  surfaces.  Integration  of  the  transformed  equations 
results  in  the  material  layer  equations. 
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Basic  Meteorological  Variables 

The  independent  variables  are  time  (t)  and  a Cartesian  system  of  space 
coordinates  (x,y,z).  The  space  coordinates  correspond  with  a position 

vector  (x1  i “ 1,2,3).  Within  this  report,  a superscripted  variable 

12  3 

denotes  a component  of  a vector.  In  this  case,  x x,  x y,  x z. 
The  coordinate  system  rotates  with  the  earth  such  that  a space  point 
remains  stationary  relative  to  a fixed  point  on  the  earth's  surface. 

The  vertical  coordinate  (x  z)  is  proportional  to  the  geopotential  (4>)* 
Thus,  z is  the  geopotential  height  with  the  origin  at  mean  sea  level. 

The  horizontal  coordinates  (x,y)  correspond  with  the  appropriate 
Universal  Transverse  Mercator  (UTM)  grid.  In  a local  region,  the  UTM 
conformal  projection  has  very  little  distortion  so  that  the  mapping 
scale  factor  is  taken  to  be  unity. 

The  principal  dependent  variables  are:  the  thermodynamic  state 
variables,  consisting  of  the  pressure  (p),  temperature  (T) , density  (p), 
and  specific  humidity  (q) ; and  the  kinematic  variables,  consisting  of 

the  velocity  vector  (u*|i  = 1,2,3;  or  u,v,w) , the  momentum  (or  mass 

flux)  vector  (m* ! i - 1,2,3,;  sometimes  written  as  mx,  my,  pT) , and  the 

if  x y 2 

water  flux  vector  (n  |i»l,2,3;  n ,n  ,n  ).  The  variables  above  are 
taken  to  be  locally  averaged  variables.  Whenever  there  is  turbulence 
present,  additional  turbulence  terms  must  be  defined.  The  more  impor- 
tant turbulence  terms  are  the  turbulent  heat  flux  vector  (r,*|i  « 1,2,3), 
the  turbulent  mass  flux  vector  (p*|i  « 1,2,3),  the  turbulent  water  flux 

vector  (v*|i  - 1,2,3),  the  Reynolds  stress  tensor  (r*^|i,j  = 1,2,3),  and 
the  turbulence  pressure  (e) . Formal  definitions  of  these  and  other 
variables  may  be  found  in  the  Glossary  of  Symbols.  Note  here  the 
following  relationships. 


i 

i i 

(la) 

m “ 

pu  + u 

i 

n - 

qm1  + v1  i « 1,2,3 

(lb) 

i 

e “ 

-l^T 

(lc) 

The  latter  equation  arises  from  the  assumption  that  turbulent  pressure 
fluctuations  are  negligible. 


Basic  Equations 

The  functional  basis  for  this  development  is  a rather  standard  system 
of  equations  for  the  state  and  kinematic  variables  in  a turbulent 
atmosphere.  The  foremost  expression  is  the  equation  of  state 

p - pRT  (2) 
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The  remaining  expressions  are  dynamic  equations.  The  simplest  is  the 
equation  of  continuity. 


■»  m.' 


o . 


(3) 


The  subscripts  in  this  and  the  following  equations  denote  partial 
differentiation.  For  example,  at  corresponds  with  ' • / • t . Whenever  re- 
peated Creek  indices  occur,  the  Klnstein  summation  convention  is  to  be 
used.  For  example. 


1 1 2 2 3 3 

m “ >m  /.>x  + 3b‘/,ix*  + am  />x 


(4) 


It  is  common  practice  for  meteorologists  to  use  the  shallow  atmosphere 
approximation  when  deriving  the  equations  of  motion  for  a rotating 
earth.  Phillips  l 20 1 has  transformed  the  equations  of  motion  from  an 
inertial  reference  frame  to  a spherical,  rotating  coordinate  system  in 
a form  which  is  consistent  with  the  shallow  atmosphere  approximation. 
(This  point  is  discussed  further  by  Veronis  [”2ll  and  Phillips  [22].) 
The  equations  of  Phillips  are  modified  in  this  report  to  conform  with 
the  UTM  Cartesian  coordinate  svstem. 


Since  viscous  effects  are  negligible  in  comparison  with  turbulence 
effects,  the  fluid  is  presumed  to  be  inviscid.  The  equations  of  motion 
govern  the  behavior  of  the  momentum: 


* + [u'u*  - t ’*  + (i  + p)iS  + m 'e  + - 0 (5) 


ii3 


\i 


where  f is  the  Coriolis  parameter;  s is  the  Kroneker  delta  which  has 


C lie  properties  that  - 1,  1 ■ 0,  for  > + i;  and  € l’  is  the  skew- 

rami 
>13 


>>i 


.'13 


123  '’13 

symmetric  tensor  which  has  the  properties  that  € ■>  1,  €~  ■ -1,  and 


t * ■ 0 for  all  other  combinations  of  Indices. 

The  moat  complicated  equation  is  that  for  the  total  energy: 


[,e  + l/2(u^u^  + 3.  )]  ^ + [me  + u'p  + 3/2u  > - i u\i  + l'  + cy  ],LC  4 Q 

(b) 


where  e is  the  sum  of  the  specific  kinetic,  geopotential,  and  internal 
energies,  Q is  the  rate  at  which  heat  is  added  per  nnit  volume  by 
radiation  and  chemical  processes,  and  I.C  is  the  rate  at  which  latent 
heat  is  added  by  condensation  of  water  vapor. 


The  continuity  equation  for  water  vapor  is  complicated  by  possible 
phase  changes  (condensation/evaporation)  and  the  precipitation  process. 
The  simplest  expression  is  that  of  the  pseudoad iabat tc  process: 


C 4 (oq)t  4 n''  - 0. 

where  C is  the  rate  ot  eondensat ion . 


17) 


In  this  simple  expression,  no  attempt  is  made  to  conserve  the  nonvapor 
phases,  which  is  equivalent  to  presuming  immediate  precipitation  ol  all 
condensed  water. 


The  basic  equations  which  have  been  enumerated  contain  numerous  esoteric 
turbulence  terms  which  have  been  included  tor  the  sake  of  completeness. 
In  practical  applications  of  these  equatiosa,  many  of  these  terms  can  be 
ignored  as  they  are  insignificant  in  comparison  with  other  terms. 

Coord lnate  Transf ormat ion 

The  geopotential  surfaces  are  not  necessarily  the  most  useful  for  del  fir- 
ing the  vertical  coordinate.  Kasahara  (23)  has  recently  reviewed  the 
general  class  of  coordinate  transformations  commonly  used  in  atmospheric 
dynamics.  Since  the  subject  of  this  report  is  the  dynamics  of  material 
layers.  It  Is  essential  that  a coordinate  system  be  chosen  so  as  to  be 
consistent  with  mass  conservation.  As  previously  noted,  there  are  natu- 
rally occurring,  quasi-horizontal,  material  surfaces  within  the  atmo- 
sphere which  may  be  used  for  defining  a nonor thogonal  coordinate  system 
which  Is  naturally  consistent  with  mass  conservation. 

Let  H(t,x,y)  and  h(t,x,y)  define  the  height  of  two  quasi-horizontal  sur- 
faces in  the  atmosphere  such  that  H h.  Between  these  surfaces  is  an 
atmospheric  layer  of  thickness  D(t,x,y)  H(t,x,y)  - h(t,x,y)  0.  For 
fixed  horizontal  coordinates,  a new  vertical  coordinate  is  defined: 

o lz  - h(t,x,y)l/D(t,x,y)  (8) 

where  h • z H.  If  the  surfaces  of  constant  sigma  are  used  as  the  ver- 
tical coordinate,  the  coordinate  system  is  no  longer  orthogonal.  Conse- 
quently, our  differential  equations  must  be  altered  to  accommodate  the 
change  of  coordinate. 

From  simple  geometrical  considerations,  the  following  rules  can  be  de- 
r ived  (see  [23]). 


CzU°i’ 


1,2 


(«a) 


or 

DC  = (DC)  I + (CDo  ) ,1-1,2  (%) 

1 1 1 v’  10 

where  C.  is  the  horizontal  derivative  of  the  arbitrary  variable  C in 
the  i-th  direction  along  a geopotential  surface,  is  the  partial 

derivative  along  the  sigma  surface,  and  D>' ^ pertains  to  the  slope  of 

the  sigma  surface  in  the  i-th  direction  relative  to  the  geopotential 
surface. 


i 


•I 

5 


10 


Do . = -(h.  + 'U.),  i - 1,2  (It 

1 1 l 

Now,  consider  the  divergence  of  any  arbitrary  dependent  vector  (a1  i = 
1,2,3).  From  equation  (9b) 


Da  | = [(Da*)x  + (Da 


%]l. 


+ an  ' 


where 


as  az  + E(axox  + ayoy)  . (] 

The  variable  as  is  equal  to  the  component  of  the  vector  normal  to  the 
sigma  surface  divided  by  cosine  psi,  where 


— X V 

P tan  D(a  o + a 
L x 


%)]  / (aXa 


X + aV)**  . 


The  change  of  vertical  coordinate  also  necessitates  a change  in  the 
definition  of  vertical  motion.  Formally 

Do  = w - (h  + oD)  , (14) 

where  the  superior  dot  denotes  the  individual  (or  substantial)  time  de- 
rivative. The  corresponding  equation  for  the  locally  averaged  "vertical" 
momentum  is 

w Doo  = m‘‘  - p(h  + oD)  = ir.S  - p(ht  + oDt)  . (15) 

This  last  equation  is  the  key  expression  of  this  report.  If  the  sigma 
surface  is  a true  material  surface,  that  is,  a surface  through  which 
there  is  no  mass  flux,  then  w(o)  = 0.  Actually,  there  is  no  assurance 
that  omega  vanishes  for  all  sigma  values,  but  u>  (1)  = 0 and  u)(0)  = 0 if 
the  layer  is  bounded  by  true  material  surfaces.  For  surfaces  which  are 
material  surfaces  in  regard  to  transport,  but  not  turbulent  diffusion, 

one  has  o = 0 and  w = PS.  Table  1 enumerates  the  basic  equations  follow- 
ing the  transformations  to  the  sigma  coordinate  system. 

Material  Layer  Integral  Operator 

Let  C denote  any  arbitrary  dependent  variable  or  expression.  The  simple 
material  layer  operator  is  defined  to  be  the  vertical  integral,  for  fixed 
(x,y),  of  the  integrand  C over  the  thickness  of  the  material  layer  be- 
tween h and  H;  that  is. 


/i’.(t,x,y)  1 

C(t,x,y,z)dz  = D(t,x,y)  J C(t,x,y,o)da  (16) 

0 
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1.  Dpt  + [ (Dmx )x  + (Dm*)  ]|o  + = 0 

2.  Dmx  + {[D(mxu  - xxx  + e + p)]  + [D(myu  - xxy)]  }| 

t x y o 

- fDmY  + [msu  - txs  + (c  + p)Do  ] =0 

' K x a 

3.  Dmjf  + {[D(mxv  - xxy)]x  + [D(myv  - Ty^  + e + p)]  }|o 

+ fDmx  + [msv  - xys  + (e  + p)Day]a  = 0 

4.  Dmz  + {[D(mxw  - xXZ)]x  + [D(myw  - xyz)]y}|o 

+ p4>a  + [msw  - xzs  + e + p]  = 0 

5.  D[pe  + l/2(u^u^  + 3e)]  + {D[mxe  + up  + 3/2ue  - xx^u^  + rx  + c cx]}  I 

t v J x1 

+ {D[mye  + vp  + 3/2vc  - xy^u^  + ry  + c cyl)  I 

v y 

+ {mse  + usp  + 3/2usc  - xs^uB  + rs  + c cs} 
v • v o 

= D(Q  + LC) 

6.  DC  + D(pq)  + [(Dnx)  + (Dny)l!  + n*  - 0 


Furthermore,  the  more  general  material  layer  moment  operator  Is  defined 
to  be 


1 

<c|n>  d/  Co "do  (17) 

0 

where  n is  a number  greater  than  zero. 

It  Is  worthwhile  to  consider  certain  special  cases.  Making  vise  of 
Leibnitz's  rule,  note  that 

<Ct>  = <C>t  - C(t,x,y,H)Ht  + C(t,x,y,h)ht  . (18) 

A similar  expression  applies  to  the  moment  operator.  Because  the  deriva- 
tives are  along  the  sigma  surfaces,  it  is  obvious  that 

<(DO1|0>  - D<C>lf  i - 1,2  (19) 

A similar  expression  applies  to  the  moment  operator.  Finally, 

<C  ,>  - D|C(t,x,y,H)  - C(t,x,y,h)  I (20a) 

<C0|n>=  DC ( t , x , y , H ) - n<c| m-1>  (20b) 

The  application  of  the  integral  operators  to  the  equations  of  table  1 
result  in  the  material  layer  equations. 

Material  Layer  Equations 

Table  2 enumerates  the  complete  elementary  material  layer  equations. 

The  normal  momentum  (w)  at  the  two  bounding  surfaces  is  included  in  tin- 
equations  for  completeness.  However,  these  terms  vanish  when  five  bound- 
ing surfaces  are  material  interfaces,  which  results  in  a significant 
simplification  of  the  equations. 

Because  of  equation  (20b),  moment  equations  always  couple  with  the  next 
lower-order  moments.  Table  3 presents  examples  of  this  coupling  as  mani- 
fested by  the  first-order  moment  material  layer  equations.  Although  the 
momentum  normal  to  the  sigma  surfaces  may  vanish  at  the  bounding  surfaces 
for  true  material  lavers,  there  is  no  a priori  assurance  that  this  vari- 
able vanishes  throughout  the  interior  of  the  layer.  Consequently,  terms 
such  as  <u'u>,  Y>,  etc.,  significantly  complicate  the  moment  mate- 

rial layer  equations. 

SOME  PRACTICAL  CONSIDERATIONS 

If  the  formal  theory  of  material  layers  were  the  only  goal,  tables  2 
and  3 would  suffice.  However,  any  attempt  to  apply  these  equations  would 
require  additional  relationships  and/or  approximations  to  effect  closure. 


13 


TABLE  2.  ELEMENTARY  MATERIAL  LAYER  EQUATIONS 


'-p'v  ♦ <mx  • + <n>y>  + ui(l)  - w(O)  = 0 

t x y 

-mx''  + vmxu  _ ixx  + t + p>  + - n/u  - ryXs  - f<my  - 

t x y 

+ [uni  - isx  + (t  + p)Dox  - Ux(ht  + oOt ) ]q  * 0 

<m^Nt  + <mxv  - ixy^x  + <myv  - iyy  + c + pN^  + f<mx> 

+ [(ov  - isy  + (e  + p)Do  - My(ht  + oDt)]Q  = 0 

-;m2>  + v.mxw  _ |XZ^  + vmYw  _ jY2''  + <p>q 

t x y J 

+ [ww  - isz  + i + p - Mz ( ht  + oD^)]^  = 0 

<pe  + l/2(i^'ul'  + 3t  )>  + <mxe  + up  + 3/2uc  - + I’x  + c <-.x' 

t v x 

+ cni^e  + vp  + 3/2vi  - i^V'  + l’y  + c fy> 

v y 

+ [uie  + usp  + 3/2usr  - ♦ Is  + CV‘*S 

- (uPuP  + 3tr ) • ( h t + 0Dt)]J,  = -Q  + LC> 

v-pqst  + >-nX>x  + <ny"  + [u>q  + vS]g  = -<CN 


♦ l .*U  - ' 


SX 


♦ (■  ♦ p)H 


- ■ XHf  ] ' 
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♦ 

[,.v  - 

isy  ♦ ( 4 P)h 

- i*yH.  ]°* 
y t 
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mxe 

4 up  + 3/2ui  - 

i xpup  + rx 

4 

• nv^e 

4 vp  + 3/2v>  - 

iy<V<  4 ry 

4 

[(.«  ♦ 

usp  ♦ 3/2us>  - 

4 rs 

4 

c ,s 

V 

- l/2(i.‘;uS  ♦ 3. 

)*(ht  + oOt 

v x 


v y 


Q ♦ LC | 1 <•  «oe  + usp  + 3/2us 


- is(V  ♦ rs  ♦ vsVD 


q|  I f nx|l  + nyH  > 0>q  + vs]'’  1 * 
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Closure  is  an  imprecise  process  which  develops  by  degrees  depending  on 
the  nature  of  the  problem  addressed.  There  are  any  number  of  specialized 
applications  for  which  the  material  layer  equations  are  appropriate,  but 
each  involves  a somewhat  different  closure  process.  This  section  con- 
siders some  of  the  practical  closure  steps  that  would  be  common  to  many 
applications. 


Turbulence  Parameterization 

In  the  context  where  they  occur,  a number  of  the  turbulence  terms  are 
negligible.  The  turbulence  pressure  is  insignificant  in  comparison  with 
the  pressure  from  the  molecular  motion.  The  horizontal  turbulent  mass 

flux  (nx,u^)  and  heat  flux  (ir,x/,^)  are  negligible  in  comparison  with  the 
corresponding  advective  fluxes.  Similar  arguments  can  be  used  to  elimi- 
nate other  turbulence  terms;  however,  the  Reynolds  stress  and  vertical 
turbulent  heat  and  moisture  flux  terms  may  be  significant  factors,  espe- 
cially for  a shallow  material  layer  in  juxtaposition  with  the  earth's 
surface. 

The  turbulent  heat  flux  (f.8)  at  the  earth's  surface  is  the  principal 
source  of  diabatic  effects  on  the  atmosphere.  The  turbulent  moisture 

flux  (vs)  is  likewise  the  principal  source  of  atmospheric  water,  and  the 

Reynolds  stress  terms  (txs,iy°)  at  the  earth's  surface  are  the  principal 
source  of  frictional  effects  on  the  atmosphere.  The  parameterization  of 
these  surface  boundary  variables  is  a closure  problem  common  to  many 
atmospheric  models.  Mel lor  and  Yamada  [24]  have  reported  on  a hierarchy 
of  turbulence  closure  models  for  planetary  boundary  layers,  and  Yamada 
and  Mellor  [25]  have  discussed  the  validation  of  one  of  these  models. 
Bhumralker  [26]  has  recently  reviewed  several  different  schemes  for 
boundary  layer  parameterization  in  atmospheric  general  circulation  models, 
and  Arya  [27]  has  further  discussed  some  of  these  approaches.  Analogous 
closure  schemes  can  be  developed  for  the  material  layer  equations,  but 
the  specific  details  are  beyond  the  scope  of  this  report. 

There  are  several  general  considerations  concerning  turbulence  effects 
on  material  layers.  Material  layers  which  are  above  the  boundary  layer 
usually  have  negligible  turbulence  effects  (free  atmosphere).  Further- 
more, it  may  be  assumed  that  turbulence  terms  are  negligible  at  the  upper 
interface  of  a material  layer  even  though  it  contains  the  boundary  layer; 
in  some  cases,  this  may  be  true  even  if  the  upper  interface  is  not  a true 
material  surface  (entrainment).  Note  that  as  a consequence  of  this 

assumption  the  principal  turbulent  dissipation  term  (ishu;  ) is  ineffec- 
tual in  the  total  energy  equation  (equation  5 of  table  2)  because  TSh 
is  assumed  to  vanish  at  the  upper  interface  and,  of  course,  u vanishes 
at  the  earth's  surface  because  of  the  no-slip  condition. 
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Although  one  may  argue  that  the  horizontal  Reynolds  stress  terms  are 
negligible  in  comparison  with  the  advective  momentum  flux,  there  are 
pragmatic  precedents  which  suggest  that  useful  approximations  of  these 
turbulence  effects  can  be  made  by  the  introduction  of  pseudoviscosity 
terms  in  the  momentum  equations.  Consistent  with  this  assumption  would 
be  the  incorporation  of  a pseudoviscosity  term  into  the  energy  equation 
to  account  for  diffusion  of  kinetic  energy. 

The  Hydrostatic  Approximation 

For  many  practical  applications  of  the  material  layer  equations,  it  is 
natural  to  resort  to  the  usual  hydrostatic  approximation  to  aid  the 
closure  process.  However,  having  eliminated  the  dynamic  equation  govern- 
ing the  vertical  momentum,  additional  approximations  must  be  introduced 
to  allow  evaluation  of  this  important  term.  The  most  forthright  approach 
to  this  problem  is  to  presume  that  all  sigma  surfaces  are  material  sur- 
faces, that  is,  that  the  momentum  normal  to  the  sigma  surfaces  (w)  is 
everywhere  negligible.  There  are  situations,  especially  in  complex  ter- 
rain, where  the  conventional  hydrostatic  approximation  and  the  latter 
premise  are  incongruous,  but  there  is  no  straightforward  substitute  for 
the  hydrostatic  approximation  which  rectifies  this  problem. 

The  assumption  of  negligible  momentum  normal  to  the  sigma  surfaces  within 
a material  layer  is  a key  consideration  in  the  closure  process  of  the 
material  layer  equations,  especially  in  regard  to  the  moment  equations. 
Table  4 enumerates  the  elementary  material  layer  equations  which  have 
been  simplified  by  the  general  closure  approximations  considered  to  this 
point.  Likewise,  table  5 contains  the  simplified  first-moment  equations. 

In  these  tables,  E 2 pe  and  K = 

Prescribed  Profiles 

The  material  layer  equations  are  nonlinear,  that  is,  they  contain  higher- 
order  terms  involving  products  of  variables.  Consequently,  some  estimate 
of  the  vertical  profile  structure  of  the  variables  must  be  available  to 
evaluate  nonlinear  terms.  The  matter  of  vertical  profiles  is  not  just  a 
pedantic  concern  for  detail.  Profile  interactions  associated  with  non- 
linear terns  affect  overall  behavior  of  the  material  layers.  In  partic- 
ular, the  density  (temperature)  profile  affects  the  buoyant  stability 
of  the  material  layer  which,  in  turn,  affects  the  kinematic  variables. 
Unlike  the  shallow-fluid  equations,  the  material  layer  equations  recog- 
nize the  compressibility  of  the  atmosphere  and  the  importance  of  thermo- 
dynamic effects.  To  capitalize  on  the  improved  generality  of  the  mate- 
rial layer  equations,  one  must  be  concerned  about  vertical  profile  struc- 
ture, even  if  only  rudimentarily. 

Assume  that  the  vertical  profile  of  any  state  or  kinematic  variable 
within  a material  layer  is  smooth  and  can  be  approximated  by  a Taylor 
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TABLE  4.  SIMPLIFIED  ELEMENTARY  MATERIAL  LAYER  EQUATIONS 


1.  <p>  + <mx>  + = 0 

t x y 

2.  <mx"t  + <mVx  + <ipyu  -y  + -p>x  - p(  1 )Hx 

+ p(0)hx  - f<my>  + ixs(0)  = kV"<.mx> 

3.  <my>^  + <mxv>x  + <myv>^  + - p(l)H 

+ p ( 0 ) + f<mx>  + ,VS (0)  = kV2<my> 

4.  <p>g  + p(l)  - p ( 0 ) = 0 

5.  <E>  ♦ p(l )H  - p ( 0 ) h + <u(E  + p)> 

t t t x 

+ <v(E  + p)>^  - <y.s(o)  = kV;,<K>  + <Q  + LC" 

6.  <pq>  + <nx>  + <„y^  _ vs(o)  = -<C> 

l x y 


TABLE  5.  SIMPLIFIED  FIRST-MOMENT  MATERIAL  LAYER  EQUATIONS 


1.  <p|l>  + <mx|i>  + <nvY|l>  =0 

t x y 

2.  <mx 1 1 >t  + <mxu|l>x  + <myu 1 1 >y  + <p 1 1 >x  - p(l)Hx  . f <my 1 1 > 

= kV:<mx | 1 > - (<p>h  + <p | 1 >D  )/D 

3.  <my|l>t  + <mxv | 1 >x  + <myv|l>y  + <p|l>y  - p(l)Hy  + f<mx|l> 

= kV-:<my 1 1 > - (<p>hy  + <p 1 1 >Dy )/D 

4.  <p| l>g  + p(l ) = <p>/D 

5.  <E|1>  + p(l)H  + <u(E  + p ) | 1 > + <v(E  + p ) | 1 > 

^ t x y 

= kV2<K| 1>  + <Q  + LC | 1 > + [c  <cS> 

P 

+ <p>h  + <p | 1 >D  - <xseu6>]/D 

v L 

6.  <pq|l>t  + <nx | 1>X  + <ny|l>y  = -<c|l>  + <vs>/D 
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series  with  "negligible"  remainder.  Consequently,  the  profile  can  be 
represented  by  a polynomial  of  some  finite  degree,  say  n.  This  approach 
parallels  that  used  for  the  development  of  multilevel  numerical  atmo- 
spheric models  where  such  Taylor  expansions  are  used  to  derive  finite- 
difference  analogs  of  the  partial  derivatives  of  the  model  equations. 

The  end  result  of  such  numerical  models  is  the  expression  of  the  verti- 
cal distribution  in  terms  of  discrete  values  at  each  level  in  the  numer- 
ical model.  The  continuous  profile  may  be  constructed  from  these  discrete 
values  by  Newton  or  Lagrange  interpolation.  On  the  other  hand,  if  one 
applies  the  material  layer  operator  for  n + 1 moments  to  the  n-degree 
polynomial,  a system  of  linear  equations  is  derived  which  relate  the 
profile  coefficients  to  the  moments.  Thus,  on  the  one  hand  the  model 
resolution  is  determined  by  the  number  of  levels  and  on  the  other  by 
the  number  of  moments  evaluated. 

Profile  description  does  not  have  to  be  restricted  to  polynomials;  any 
integrable  function  can  be  used.  For  example,  since  a power-law  is 
commonly  used  to  approximate  the  wind  profile  in  the  boundary  layer, 
this  would  be  a more  logical  choice  than  a zero-  or  first-order  polynomial. 
On  the  other  hand,  the  combination  of  the  hydrostatic  approximation  and 
the  equation  of  state  restricts  the  degrees  of  freedom  for  character- 
izing the  profiles  of  state  variables.  The  simplest,  nontrivial  case 
is  a polytropic  representation. 

To  this  point,  only  profiles  within  a material  layer  have  been  discussed. 

If  the  atmosphere  is  composed  of  several  overlapping  material  layers,  the 
complete  profile  is  pieced  together  from  the  profiles  of  each  layer.  At 
the  interface  between  material  layers,  the  dynamic  boundary  condition 
must  be  satisfied.  (The  kinematic  boundary  condition  is  met  a priori.) 
Indeed,  it  Is  natural  to  require  continuity  of  all  variables  across  the 
interface.  However,  there  are  practical  reasons  for  relaxing  the  require- 
ment of  continuity. 

Consider  a hypothetical  case  where  the  atmospheric  structure  in  a given 
region  consists  of  a single  subsidence  inversion  imbedded  in  an  other- 
wise comparatively  labile  troposphere.  Most  likely,  the  vertical  dis- 
tribution of  moisture  and  pollution  would  reveal  that  there  is  a material 
layer  below  the  subsidence  inversion.  Likewise,  the  layer  between  the 
subsidence  inversion  and  the  tropopause  would  be  another  material  layer. 
Profiles  within  these  layers  could  be  approximated  by  low-order  poly- 
nomials provided  no  portion  of  the  subsidence  inversion  is  incorporated 
into  the  layers.  To  complete  the  continuous  profile,  it  is  necessary  to 
include  the  subsidence  inversion  layer  as  a third  material  layer  buffer 
between  the  upper  and  lower  layers. 

Numerous  situations  arise  as  described  above,  where  comparatively  thick 
material  layers  are  separated  by  a stable  layer  which  is  thin  but  still 
significant  in  terms  of  profile  structure.  Because  of  the  almost  insig- 


nif leant  thickness  of  these  layers,  it  is  desirable  to  approximate  these 
layers  by  an  infinitesimally  thin  interfacial  layer.  In  considering 
the  dynamic  boundary  condition  in  relation  to  the  equations  of  table  4, 
the  following  diagnostic  equation  has  been  derived  to  govern  the  inter- 
facial values. 

|l/2(n  u.'.u  + ir.XvAv)  + (c^/R)pNu  + c ih*AtJ  + |\/2(m'uAu  + m^vAv) 


+ (Cp/R) pAv  + cpmVAT  = 0 (21) 


In  this  equation,  AT,  Au,  and  Av  are  the  magnitude  of  the  jump  discon- 
tinuities of  temperature  and  velocity  across  the  interface  and  the  supe- 
rior bar  denotes  "average"  values  of  the  variables  within  the  infinitesi- 
mal inter  facial  layer.  For  practical  purposes,  these  quantities  are  the 
average  of  the  appropriate  upper  and  lower  interfacial  values. 

The  introduction  of  interfacial  discontinuities  is  an  artifice  but  is  of 
practical  value  in  simplifying  the  modeling  of  material  layer  behavior. 
Equation  (21),  which  is  a generalized  representation  of  Margules’  equa- 
tion, allows  for  a physically  meaningful  control  over  the  magnitude  of 
the  jumps. 

The  piecewise  construction  of  vertical  profiles  generally  requires  a 
knowledge  of  the  variables,  or  quantitative  relationships  which  govern 
their  behavior,  at  the  upper  and  lower  bounding  interfaces  as  well  ns 
the  internal  interfaces. 


Bounding  Passive  Layer 

For  most  practical  applications,  it  is  neither  necessary  nor  desirable 
to  apply  the  material  layer  equations  to  the  entire  depth  of  atmospheric 
mass.  Most  problems  of  interest  are  confined  to  the  troposphere.  Under 
these  circumstances,  it  is  necessary  to  introduce  a superior  bounding 
layer  which  is  more  or  less  passive  in  relation  to  the  material  layer 
effects  which  occur  below.  Since  the  objective  is  to  establish  the  req- 
uisite boundary  conditions  at  the  uppermost  material  interface,  the  pre- 
scription of  the  structure  of  the  passive  layer  need  not  be  in  any  greater 
detail  than  that  required  to  achieve  the  objective. 

For  the  first  example  of  a passive  bounding  layer,  consider  a hypothetical 
stratosphere  which  is  isothermal.  According  to  the  hydrostatic  equation, 
the  pressure  and  density  profiles  have  exponential  decay  with  height  and 
there  is  no  distinct  "top"  for  the  atmosphere.  Since  no  thermal  wind 
exists  in  this  hypothetical  stratosphere,  it  is  reasonable  to  assume  that 
the  wind  is  constant  with  height.  It  is  readily  shown  that  knowledge  of 
the  state  and  kinematic  variables  on  the  tropopause  is  sufficient  to 
allow  specification  of  all  variables  throughout  the  stratosphere. 


Table  6 enumerates  a system  of  equations  governing  the  state  and  kinematic 
variables  on  the  tropopause  below  an  isothermal  stratosphere  with  pre- 
scribed geostrophic  wind  (u  , v ).  Here,  H and  p are  the  height  and 

P o 

pressure,  respectively,  of  the  tropopause  and  all  space  derivatives 
pertain  to  the  tropopause  surface.  The  solution  of  these  equations 
provides  the  requisite  interfacial  values  to  couple  with  the  tropospheric 
profiles . 

Further  simplifications  of  the  passive  isothermal  layer  may  be  effected 
by  assuming  the  stratospheric  wind  (including  the  tropopause)  to  be 
geostrophic.  In  either  case,  the  isothermal  assumption  may  be  crude  but 
not  unreasonable,  and  should  serve  well  for  establishing  tropopause  values 
for  a tropospheric  material  layer  model.  A more  general  model  is  achieved 
by  using  the  isothermal  stratospheric  model  to  specify  upper  interfacial 
values  in  conjunction  with  the  incorporation  of  equation  (21)  into  the 
model . 

The  simplest  possible  example  of  a passive  bounding  layer  within  the 
troposphere  is  one  assumed  to  be  adiabatic  with  the  wind  uniformly 
geostrophic.  The  bounding  interface  is  presumed  to  be  an  isentropic 
surface  and  consequently  the  Montgomery  streamf unct ion  (pH  + c(iT)  is  a 

linear  function  of  position  on  the  interface  (as  a consequence  of  the 
geostrophic  wind  assumption).  The  wind  and  temperature  profiles  within 
the  interior  material  layer  must  be  defined  so  as  to  approach,  as  o * 1, 
the  geostrophic  wind  and  the  temperature  consistent  with  the  Montgomery 
streamf unction  relationship  alluded  to  above.  Again,  a more  general 
model  is  achieved  by  using  the  adiabatic  layer  model  to  prescribe  the 
upper  interfacial  values  for  use  in  equation  (21). 

The  examples  of  upper  boundary  conditions  described  above  are  only  two 
of  many  possible  alternatives;  the  choice  depends  upon  the  application 
to  be  made  of  the  material  layer  equations.  A major  point  to  note  is 
how  additional  degrees  of  freedom  are  allowed  for  the  material  layer 
solutions  by  use  of  the  generalized  Nargules  equation  (21),  which  is  a 
diagnostic  equation  and  thus  readily  used  in  numerical  models. 

Air/Earth  Interface 

Air/earth  interaction  at  the  ground  surface  is  intimately  related  to  the 
problem  of  boundary  layer  turbulence;  however,  it  is  beyond  the  scope  of 
this  report  to  present  detailed  air/earth  interaction  relationships  for 
application  to  material  layer  models.  Regardless  of  the  detailed  approach 
used,  the  kinematic  variables  must  vanish  at  the  air/earth  interface  in 
accordance  with  the  no-slip  condition,  and  the  temperature  (and  specific 
humidity)  must  be  prescribed  or  evaluated  from  an  appropriate  surface 
energy  budget  model. 

The  boundary  conditions  at  the  air/earth  interface  are  a problem  common 
to  many  atmospheric  prediction  models.  For  example,  bhumralker  [28]  has 


•)  1 


r 


TABLE  6.  TROPOPAUSE  EQUATIONS  FOR  ISOTHERMAL  STRATOSPHERE 
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described  a hierarchy  of  numerical  algorithms  for  the  computation  of 
ground  surface  temperature  in  an  atmospheric  general  circulation  model, 
and  Deardorff  [29]  has  further  developed  this  approach  to  include  a layer 
of  vegetation.  In  arriving  at  an  appropriate  air/earth  interaction  model, 
a critical  decision  must  be  made  as  to  whether  the  evaluation  of  the  soil 
heat  flux  requires  that  the  "soil  temperature"  be  prescribed.  A prescribed 
soil  temperature  provides  negative  feedback  to  an  energy  budget  model 
which  assures  reasonably  predicted  surface  conditions  (provided  the  speci- 
fied soil  temperature  is  reasonable). 

The  Bhurcralker /Deardorf f approach  can  be  generalized  for  use  with  material 
layer  equation  models,  but  for  complex  terrain  the  specification  of  the 
many  esoteric  parameters,  least  of  all  "soil  temperature,"  may  be  by  guess. 
Under  these  circumstances,  it  might  be  just  as  well  to  prescribe  outright 
the  interface  temperature  and  thus  forego  the  complexities  of  a surface 
energy  budget  model  of  questionable  validity. 

MATERIAL  LAYER  MODELS 

The  central  theme  of  this  report  is  the  theoretical  development  of  the 
material  layer  equations  and  the  associated  problems  of  closure.  Exam- 
ples of  material  layer  models  are  presented  in  this  section  to  demon- 
strate various  specialized  applications  of  the  material  layer  equations. 
Although  numerical  algorithms  are  not  presented  in  support  of  the  exam- 
ples, the  utility  of  the  sample  material  layer  models  should  be  evident. 

The  first  example  is  presented  as  the  material  layer  analog  of  the  con- 
ventional one-layer  shallow-fluid  equation  terrain  effects  model.  The 
second  example  is  an  extremely  simple  model  to  demonstrate  the  utility 
of  the  material  layer  equations  when  the  bounding  interface  is  not  a 
material  surface.  Finally,  the  third  example  deals  with  the  use  of  the 
material  layer  equations  for  diagnostic  analysis  of  terrain  effects  on 
atmospheric  phenomena  using  a variational  approach. 

An  Elementary  Mixing  Layer  Model 

Tingle  and  Bjorklund  [5]  have  described  a mixing  layer  model  based  upon 
the  shallow-f luid  equations;  the  objective  here  is  to  present  the  basis 
for  a mixing  layer  model  derived  from  the  material  layer  equations.  Al- 
though intended  as  the  compressible  fluid  analog  of  the  shallow-fluid 
equations,  a comparison  of  the  two  approaches  is  like  that  of  apples  to 
oranges.  However,  it  is  appropriate  to  consider  this  as  the  basis  for 
further  generalization  of  the  Tingle-Bjorklund  model  to  incorporate, 
among  other  things,  thermodynamic  effects  on  mixing  layer  dynamics. 

The  mixing  layer  is  defined  to  be  the  primary  material  layer  in  juxta- 
position with  the  ground.  It  may  be  a rather  deep  labile  layer  capped 
by  an  elevated  subsidence  inversion  or  a shallow  stable  layer  confined 
to  a nocturnal  surface  inversion.  Since  only  a one-layer  model  is  being 
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considered,  the  isentropic,  geostrophic  material  interface  described  in 
"Bounding  Passive  Layer"  is  chosen  as  the  upper  boundary  condition.  Ihe 
layer  structure  is  assumed  to  be  polytropic  (linear  temperature  profile). 
The  "ground"  temperature  could  be  implicitly  evaluated  if  there  was  pro- 
vision for  an  air/earth  interaction  model.  However,  for  this  elementary 
mixing  layer  model,  it  shall  be  assumed  that  the  "ground"  temperature  is 
steady  and  everywhere  prescribed  and  that  the  turbulent  heat  flux  at  the 

ground  [ Cs (0) ] is  negligible.  Table  7 enumerates  the  relationships  gov- 
erning the  layer  structure.  Table  8 lists  the  boundary  conditions  at  the 
air/earth  interface,  and  table  9 itemizes  the  boundary  conditions  at  tl 
upper  isentropic,  geostrophic  surface.  Finally,  table  10  enumerates  the 
material  layer  equations  for  the  elementary  mixing  layer  model. 

Aside  from  the  coriolis  force,  surface  drag,  and  pseudoviscosity  terms 
which  the  Tingle-Bjorklund  model  neglects,  the  first  three  equations  of 
table  10  are  analogous  to  those  of  Tingle-Bjorklund  provided  that  the 
density  is  assumed  constant  and  the  vertical  profile  of  velocity  is 

uniform.  The  term  g0(’>  1'  - 1/2,  D)  is  analogous  to  the  Tingle-Bjorklund 

term  l/2g‘D“,  which  is  the  square  of  the  gravity  wave  speed,  where 

~ o)/o  is  the  reduced  gravity.  The  Tingle-Bjorklund  model  has  no 
counterpart  to  equation  4 of  table  10. 

Several  points  should  be  made  concerning  the  derivation  of  equation  4 of 
table  10.  First,  there  are  no  time  derivatives  of  the  temperature  because 
of  the  manner  in  which  the  upper  and  lower  boundary  conditions  were  stated. 
Furthermore,  the  first  moment  of  the  continuity  equation  was  used  in 
deriving  equation  4.  However,  the  main  point  is  that,  by  comparison  with 
the  shallow-fluid  equations,  the  fourth  equation  introduces  a whole  new 
perspective  to  terrain  and  thermodynamic  effects  on  mixing  layer  dynamics. 

The  use  of  the  Isentropic,  geostrophic  upper  material  interface  would  be 
most  appropriate  for  the  case  of  a mixing  layer  associated  with  a surface- 
based  inversion.  On  the  other  hand,  if  the  upper  material  interface  is 
the  base  of  an  elevated  inversion,  a more  general  model  is  achieved  by 
the  use  of  equation  (21)  wherein  the  isentropic,  geostrophic  assumption 
pertains  only  to  the  upper  interface.  Furthermore,  equation  (21)  would 
be  modified  for  the  simplest  model  by  assuming  uniformity  in  the  inter- 
facial layer,  that  is 

l/2(mxuAu  + mxvAv)  + (cp/R)pAu  + c mxAT  = constant  (22a) 

l/2(m^uAu  + m^vAv)  + (c  /R)pAv  + c m^'.T  = constant  (22b) 

P P 

The  profile  relationships  of  table  7 must  be  modified  to  reflect  values 
at  the  upper  boundary  of  the  mixing  layer  consistent  with  the  lower 
Interfacial  values  derived  from  equations  (22a)  and  (22b).  Although 
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TABLE  7.  LAYER  STRUCTURE  FOR  ELEMENTARY  MIXING-LAYER  MODEL 


1 . Lapse  Rate 

Y = -Tz  = [T (0 ) - T(1  )]/D 

2.  Polytropic  Pressure  Profile 

p(o)  = p(0)[T(o)/T(0)]<5 , where  6 = g/Ry 

• ‘ 

3.  Linear  Temperature  Profile 

T(o)  = T(0)  - yDo 

f 

4.  Pressure  Integral  L* 

<P"  = D[g<p|l>  + p(l )]  = R[(<p>  - <p 1 1 > )T ( 0)  + <p 1 1 >T ( 1 ) ] t 

5.  Prescribed  Power-Law  Mixing-Layer  Wind  Profile 

u(o)  = u on(l-a)  + Oo 
o 

v(o)  = v on( 1-a ) + tfo 
o 

where  n is  prescribed.  ji 


TABLE  8. 


BOUNDARY  CONDITIONS  AT  THE  AIR/EARTH  INTERFACE 


1.  Temperature  Prescribed 

T(O)  = T(x,y) 

2.  No-Slip  Condition 

u(0)  = v (0)  = 0 

3.  Form  Drag  Relationship  For  Surface  Reynolds  Stress 

iXS(0)  = du  , Tys(0)  = dv  , 
o o 

where  d is  also  a function  of  uq,  vq  and  the  surface  roughness. 

4.  Negligible  Heat  Exchange 

f.S(0)  = 0 


II 


i 


- 
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TABLE  9.  BOUNDARY  CONDITIONS  AT  UPPER  ISENTROPIC  GEOSTROPHIC  SURFACE 


1.  Montgomery  Streanifunction  Relationship 

gH  + c T(i)  = (gH  + cpT)  + f(vx  - uy)  , 

where  (gfl  + c^T)  is  a nominal  (constant)  Montgomery  streanifunction, 
and  0,  are  the  (constant)  geostrophic  wind  components. 

2.  Pressure  Relationship 

P 11(1)  = P [p(l)/p  ]R/Cp  = RrT(l) 

0 0 0 

where  r is  the  potential  density  of  the  isentropic  surface. 

3.  Wind  Relationship 

u(l)  = Q , v ( 1 ) = 0 


'' 


J 


• 4 


1 
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TABLE  10.  MATERIAL  LAYER  EQUATIONS  FOR  ELEMENTARY  MIXING-LAYER  MODEL 


1.  <p>.  + + --my>  = 0 

t x y 

2.  <mx>t  + ‘-nixu''x  + + g[D(<p|l'*  - 1 /2f?D ) ]^ 

+ g[<P'>  - pD]h  - f ( 11  - PDS/)  + du 
x o 

3.  <my  * + '-mxv'>  + ni-yv>  + g[D(<p|l>  - l/2pD)] 

t x y y 

+ g[<p>  - f'D]h  + f(^m*N  - ^Du)  + dv 

y o 


kV-\mx> 


kV;>^m^N 


4.  [•  K ' - gD(<p|l>  - 1 /2pD) ] + vUK>  + <vK" 

t x y 

+ (<mx>  - ^mx 1 1 ^ ) • (gh  + cpT)x 
+ (<mys  - <my|l>)*(gh  + cpT)y 
+ f[vmx|l>y'  - Miiy  1 1 '•Q]  = kV2<K> 
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this  further  complicates  the  solution  of  the  model  equations,  additional 
degrees  of  freedom  are  provided  for  profile  evaluations  without  the  intro 
duction  of  additional  dynamic  differential  equations. 

Mixing  Layer  Entrainment  Model 

The  presence  of  an  inversion  does  not  necessarily  assure  the  existence 
of  a material  layer.  The  processes  of  penetrative  convection  and  entrain 
ment  invalidate  the  premise  of  mass  conservation.  In  the  case  of  penetra 
tive  convection,  organized  convective  elements  have  sufficient  buoyancy 
to  "punch  holes"  in  the  inversion.  The  inversion  is  more  like  a sieve 
than  a lid,  and  consequently  it  is  difficult  to  reconcile  the  material 
layer  equations  with  the  condition  of  penetrative  convection.  In  the 
case  of  entrainment,  one  or  both  layers  adjacent  to  the  inversion  are 
turbulent  with  the  result  that  the  turbulence  erodes  away  the  inversion 
layer  with  an  associated  entrainment  of  mass  at  the  expense  of  the  inver- 
sion layer. 

The  classical  example  of  entrainment  is  the  burning-off  of  a nocturnal 
surface  inversion  after  sunrise.  The  solar  heating  leads  to  a mixing 


layer  which  deepens  with  time  and  eventually  wipes 
layer.  To  illustrate  this,  consider  a simple  case 
horizontal  homogeneity.  The  behavior  of  the  mixing 
by  the  following  equations: 

out  the  inversion 
with  calm  wind  and 
l layer  is  described 

<P>t  + w(l)  = 0 

(23) 

<P>t  + RMl)T(l)  + ts(l)]  = Rts(0) 

(24) 

where  -u>(l)  is  the  rate  of  mass  entrainment,  4s  (1)  and  4s  (0)  pertain  to 
the  turbulent  heat  flux  at  the  bottom  and  the  top  of  the  mixing  layer, 
respectively.  If  the  layer  above  the  inversion  is  a nonturbulent  mate- 
rial layer,  the  combined  mixing  layer  and  inversion  layer  is  also  a 
material  layer  (this  is  the  difference  between  entrainment  and  penetra- 
tive convection) . The  material  layer  equations  for  the  combined  mixing/ 
inversion  layer  are 

<P>t  « 0 

(25) 

<P>t  - Rf S(0) 

(26) 

This  formulation  of  the  problem  should  be  compared  to  that  given  by 
Tennekes  [30].  To  facilitate  this  comparison,  it  is  assumed  that  the 
mixing  layer  has  a dry  adiabatic  lapse  rate.  Consequently 


<e> 


Mixing  layer 


(P(O)T(O)  - p(l) . T (1) ) / (g/R  + g/c  ). 


(27) 


The  interior  of  the  inversion  layer  is  assumed  nonturbulent  and  adiabatic 
with  respect  to  time  (up  to  the  time  the  inversion  is  annihilated).  For 
the  sake  of  simplicity,  assume  that  the  lapse  rate  within  the  inversion 
layer  is  a constant.  Because  of  the  hydrostatic  and  adiabatic  assumptions, 
the  lapse  rate  within  the  inversion  layer  is  constant  in  time  as  well  as 
space.  It  follows  from  the  preceding  that 

rS(l)  - 0 (28) 


and 


<*>(1)  - P(1)Cb(0)/(y  - 


R/c  ) 

P 


<P>, 


Mixin",  l aver 


(29a) 


« rS(0)/(>  - g/c  )D  (29b) 

P 

where  Y is  the  lapse  rate  in  the  inversion  layer.  Combining  equations 
(23)  and  (29b),  one  may  derive  the  approximate  relationship 

(D2)  a -2CS(0)/P(>  - g/c  ) (30) 

t P 

where  >’  is  the  mean  mixing  layer  density.  For  a steady  surface  heat  flux, 
equation  (30)  is  in  consonance  with  the  square-root-of-time  relationship 
for  the  evaluation  of  the  mixing  laver  depth  (see  equation  (19)  of  Tennekes 
[30]). 

The  material  layer  approach  to  the  problem  of  daytime  mixing  layer  evolu- 
tion is  quite  formal  but  straightforward.  There  are  significant  differ- 
ences between  this  approach  and  that  of  Tennekes  [30],  but  comparison  is 
difficult.  Ter.nekes  speaks  of  entrainment  and  penetrative  convection  as 
if  they  were  the  same  process,  possibly  because  he  introduces  the  artifice 
of  characterizing  the  inversion  layer  as  a jump  discontinuity  rather  than 
a layer  of  finite  thickness. 

The  elementary  model  of  daytime  nixing  layer  evolution  given  above  is  one 
example  of  the  utility  of  the  material,  layer  equations  even  when  the  con- 
dition of  mass  continuity  is  negated.  Similar  "leaky"  material  layer 
applications  are:  the  opposite  problem  of  the  evolution  of  the  nocturnal 
surface  inversion  layer,  the  more  general  problem  of  nocturnal  air  drain- 
age, and  the  related  problem  of  diurnal  slope  winds  associated  with  sur- 
face heating. 


Variational  Material  I. aver  Model 

Once  suitable  lateral  boundary  conditions  are  established,  the  mixing 
layer  model  set  forth  in  "An  Elementary  Mixing  Layer  Model"  becomes  a 


deterministic  initial-value  problem  which  is  most  suitable  for  prognos- 
tication. Even  so,  models  of  this  type  have  also  been  used  for  diagnostic 
analysis  [3,4,5].  However,  there  are  definite  problems  in  the  use  of 
wholly  deterministic  models  for  diagnostic  analysis  since  any  attempt  to 
constrain  the  model  solution  to  conform  to  observational  input  data 
(except  for  boundary  and  initial  conditions)  may  lead  to  an  ill-posed 
problem.  In  this  case,  variational  analysis  may  be  used  to  draw  a com- 
promise between  model  equations  and  observational  input  data  for  diag- 
nostic analysis. 

Ohmstede  [31]  and  Meyers,  Cedarwall,  and  Ohmstede  [32]  have  reported  on 
a model  for  diagnostic  analysis  of  terrain  effects  on  atmospheric  phe- 
nomena using  a variational  approach.  This  model  is  so  straightforward 
that  its  principles  are  readily  summarized  here.  At  any  one  time,  the 
model  deals  with  a single  material  layer  of  thickness  D(x,y)  > 0 in 
juxtaposition  with  the  ground.  The  material  layer  variables  subjected 
to  variational  adjustment  are  the  depth  of  the  layer,  the  zero-  and 
first-order  moments  of  the  mass  flux  vector,  the  total  energy  flux  vec- 
tor, and  the  potential  energy  flux  vector.  The  height  of  the  ground  is 
prescribed. 

The  general  concepts  of  the  variational  analysis  begin  with  some  simply 
connected  region  for  which  analysis  is  desired.  Using  the  limited  obser- 
vations available,  some  form  of  objective  analysis  is  used  to  construct 
fields  of  the  relevant  variables  within  the  region.  These  are  called 
the  "observed"  fields  even  though  the  only  observations  are  actually  at 
discrete  space/tirae  points.  Only  in  a very  limited  way  do  the  observed 
fields  reflect  terrain  effects.  Consequently,  the  desire  is  to  construct, 
by  means  of  variational  analysis,  new  fields  which  better  characterize 
the  terrain  effects  but  still  are  in  consonance  with  the  observed  fields 
in  a least-squared  sense.  The  object  is  to  minimize  the  square-variation 
of  the  field  variables  but  subject  to  functional  constraints  which  mani- 
fest terrain  effects.  The  constraints  used  are  the  dynamic  equations  1 
and  5 of  table  4 and  equation  1 of  table  5,  and  two  algebraic  equations 
which  characterize  the  interdependence  of  the  variables.  The  resultant 
variational  fields  of  this  analysis  model  are  the  observed  fields  adjusted 
to  comply  with  the  material  layer  constraints. 

Many  different  variational  models  are  possible  using  the  material  layer 
equations.  More  advanced  models  should  consider  multilayer  interactions 
and  the  use  of  the  momentum  equations. 

CONCLUSIONS 

The  objective  of  this  report  has  been  to  set  forth  new  tools  for  atmo- 
spheric diagnostic/prognostic  analysis.  It  could  be  argued  that  the 
material  layer  equations  are  not  really  new  but  just  a further  general- 
ization of  the  shallow-fluid  equations.  Although  largely  true,  the 
material  layer  equations  contain  unique  features  which  are  entirely  over- 
looked by  the  incompressible  shallow-fluid  approach. 
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The  material  layer  approach  is  most  useful  when  the  matter  of  concern  is 

the  transport  and  diffusion  of  material  such  as  from  NBC  weapons  or  smoke. 

Equally  important,  however,  are  the  features  of  the  material  layer 

approach  which  allow  a more  general  analysis  of  the  interaction  of  ter- 

rain and  thermodynamic  effects.  Although  the  material  layer  equations 
could  surely  be  used  for  prognosis  if  suitable  initial  and  boundary  con- 
ditions could  be  established,  it  is  believed  that  their  greatest  utility 
is  likely  to  be  for  mesoscale  diagnostic  analysis  of  terrain  and  thermo- 
dynamic effects  using  a variational  approach. 
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Watts/m' 

kgm/m' 

j/kgm-deg 


sec 
deg  K 
deg  K 


Drag  coefficient 

Specific  heat  at  constant  pressure 

Specific  heat  at  constant  volume 

Surface  stress  coefficient 

Depth  of  material  layer 

Specific  total  energy 

Volumetric  total  energy 

Coriolis  parameter 

Acceleration  of  gravity 

Elevation  of  base  of  material  layer 

Elevation  of  top  of  material  layer 

Pseudoviscosity 

Kinetic  energy 

Heat  of  vaporization 

Momentum  component  in  i-th  direction 

Water  vapor  flux  component  in  the 
i-th  direction 

Pressure 

Reference  pressure 

Specific  humidity 

Heating  rate  due  to  radiation  and 
chemical  processes 

Potential  density 

Gas  constant 

Superscript  denoting  coordinate 
direction  normal  to  u-surfaces 

Time 

Temperature 
Ground  temperature 
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Turbulent  velocity  fluctuation  in 
i-th  direction 

Velocity  component  in  X-direction 
Velocity  component  in  Y-direction 
Velocity  component  in  Z-direction 
Coordinate  in  i-th  direction 
UTM  grid  easting 
UTM  grid  northing 
Elevation 

Coordinate  index  requiring  summation 
convention 

Coordinate  index  requiring  summation 
convention 

Lapse  rate 

Turbulent  transfer  of  turbulent 
kinetic  energy 

Kroneker  delta 

Turbulence  pressure 

Skew-symmetric  tension 

Turbulent  heat  flux  in  i-th  direction 

Turbulent  mass  flux  in  i-th  direction 

Turbulent  water  vapor  flux  in  i-th 
di rection 

Adiabatic  reduced  pressure 
Densi ty 

Material  layer  vertical  coordinate 
Reynolds’  stress  tensor 
Geo potential 

Momentum  component  normal  to 
o-surfaces 
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